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A novel derivation of the Q? dependence of quark and gluon densities (of given 
helicity) as predicted by quantum chromodynamics is presented. The main body of 
predictions of the theory for deep-inleastic scattering on either unpolarized or polarized 
targets is re-obtained by a method which only makes use of the simplest tree diagrams 
and is entirely phrased in parton language with no reference to the conventional operator 
formalism. 


1. Introduction 


The quark parton model [1] provides us with a very useful and simple description 
of the physics of deep inelastic phenomena [2]. The theoretical framework which 
justifies the parton model is given by the asymptotically free gauge theory of strong 
interactions based on the color degrees of freedom [3] (quantum chromodynamics, 
QCD). Although scaling is predicted to be broken by logarithms (a fact which appears 
to be well consistent with present experiments), the deviations from scaling can be 
and have been computed for deep inelastic structure functions for either unpolarized 
[4,5] or polarized targets [6,7]. In the leading logarithmic approximation, the results 
can again be phrased in the parton language by assigning a well determined Q? depend- 
ence to the parton densities. In spite of the relative simplicity of the final results, 
their derivation, although theoretically rigorous, is somewhat abstract and formal, 
being formulated in the language of renormalization group equations for the coeffi- 
cient functions of the local operators which appear in the light cone expansion for 
the product of two currents. 
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In this paper we show that an alternative derivation of all results of current inter- 
est for the Q? behaviour of deep inelastic structure functions is possible. In this ap- 
proach all stages of the calculation refer to parton concepts and offer a very illumin- 
ating physical interpretation of the scaling violations. In our opinion the present 
approach, although less general, is remarkably simpler than the usual one since all 
relevant results can be derived in a direct way from the basic vertices of QCD, with 
no loop calculations being involved (the only exception is the lowest order expres- 
sion for the running coupling constant which we do not rederive). 

This method can be described as an appropriate generalization of the equivalent 
photon approximation in quantum electrodynamics [8]. A preliminary and less 
complete version of this paper has been presented by one of us at the 1976 Flaine 
meeting [9]. 

The present paper is organized as follows. In sect. 2 the basic formulae for deep 
inelastic processes are recalled and notations and definitions are specified. In sect. 

3 the integro-differential equations that describe the Q? dependence of parton 
densities are introduced and their physical interpretation is discussed at length. In 
sect. 4 the calculation of logarithmic exponents for the spin averaged structure 
functions is described, while sect. 5 is devoted to the calculation of the other loga- 
rithmic exponents which are relevant to spin dependent structure functions. Sect. 6 
contains some concluding remarks. 


2. Basic formulae and definitions 


In this section we recall some basic notions and notations concerning deep inelas- 
tic scattering on a polarized target which are of relevance for the following. For 
simplicity in all of this paper we refer to the case of electroproduction, the exten- 
sion to neutrino scattering being simple and well known. 

Let p and q be the four-momenta of the nucleon (of mass M) and of the virtual 
photon respectively. As usual we define: 


‘Mv = (pq) , -g?=@?>0, 2Mrx=Q?. (1) 


The structure functions of electroproduction are defined as [1,6] 


E ; 
4r? —| dye! (p, s\[J 5), JAS) Ip, s) 
M 


dv Mi Mv \ W(x, Q 
= (2w k F ) W(x, Q?) + (>, -7 au) (Po a a) a 
2 2 
E ls Heg) +(e" = (ap?) aes i, é 


where s° is the polarization four-vector of the nucleon. The number of independent 
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structure functions is fixed by the number of non-vanishing independent helicity 
amplitudes for forward virtual Compton scattering in the case of parity and time 
reversal invariance. 

Standard arguments based on free quark light cone expansion [10] or the parton 
model [1] show that 


MW, =F, ? vW =F), (3) 


vV =G, vV, =MG; , (4) 


are finite and independent of Q? in the scaling approximation (that is when loga- 
rithmic effects from strong interactions are neglected). In this same approximation 
it is also true that 


2xF =F, (5) 


as a consequence of the vanishing of the longitudinal cross section. Moreover, the 
structure functions F, F2, and G, can be simply expressed in terms of parton 
quark densities. Let g‘(x) be the number density of quarks of type i (summed over 
colors) inside a proton target with fraction x of the proton longitudinal momen- 
tum in the pæ frame. Similarly we denote, by q‘ (x) and qê (x) the densities of 
parton quarks with positive and negative helicities in a proton of positive helicity. 
One has, of course, 


qi.(x) + q(x) =q"(x). (6) 
Defines e; as the charge of the ith quark in units of the proton charge, we have [1] 


2F; = Fy/x = 27 e? [i +T) , (7) 


2G, = 2 e? li -iO T LTE - (8) 


The structure function G, does not have an equally simple expression in terms 
of quark densities of given helicity. This is related to the fact that G would be 
zero for a point-like proton (namely a proton that is itself a parton), as one can 
easily see from the Born diagram of virtual Compton scattering. In light cone 
language one finds that G2 would be zero if matrix elements of local operators 
were extracted from free quark field theory. This is certainly not admissible in 
general. However our experience (for example in the case of ‘sea’ densities) strong- 
ly suggests that Gz will turn out to be a small effect. 

We recall that G, contributes [10], quite in general, corrections of order 1/Q? 
to the asymmetry in the case of longitudinally (with respect to the beam) polarized 
protons. In this case the dominant contribution of order one is determined by G4. 
For transversely polarized protons, G; and Gz contribute [10] terms of the same 
order to the asymmetry which however is itself of order 1/Q in this case. 
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In this paper we shall restrict ourselves to a discussion of scale breaking effects 
concerning F', Fa and G, which have a simplest description in parton language. 
In fact in the following sections we shall show that the Q? dependence of the densi- 
ties qi, +q! and qi‘, —q'_ for all types of quarks and antiquarks can be directly 
evaluated in QCD on the basis of simple physical arguments. From the previous 
discussion we see that this knowledge, although not complete is however sufficient 
to reproduce all the results of the theory which are of current interest. 


3. The master equations 


In this section we shall introduce a set of integro-differential equations for the 
Q? dependence of quark and gluon densities which are the basis for the present 
approach. These equations were originally obtained [11] by Mellin transform tech- 
niques starting from the renormalization group equations. Hereafter briefly recalling 
their derivation we discuss their physical contents [12]. We shall deal first with the 
unpolarized case when only the densities 


E) = ai) + gL () 
are relevant. The further extension to the general case will be done in sect. 5. 
We start by considering the simplest case of only one flavour of quarks. We de- 


note by qS the “net” number of quarks in the proton as seen by the current in the 
Po frame, that is the number of quarks minus antiquarks: 


aS(x, t) = ax, t) = qx, t) 4 (9) 
The variable ¢ is defined as 
t=In Q?/Q5 , (10) 


with Qe a suitable normalization point. The label NS stands for ‘“‘non singlet” since, 
in a more general situation, differences of this type transform as the adjoint represen- 
tation of the flavour group. The Q? dependence implied by QCD can be simply 
expressed in terms of moments of parton densities. We therefore define * 


1 
MNS) = fax x"! Se, 2). (11) 
0 


As is well known [4,5] the predicted t dependence of moments is of the form 
per 
“Nn 


MBS(1) = MNS(O) Es 


FO) om 


* In the following unless explicitly stated we restrict to moments with n > 1 in order to have 
convergence at x = 0 in moments integrals. 
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Here a(t) is the running coupling constant of QCD (a= g?/4m in terms of the 
quark gluon vertex). In the leading logarithmic approximation a(t) is of the form 
[4,5] 


=1+bat, 13 
a 
with a = o{0) and b a constant which is given in general by 

zi C,(G) — 4 T(R) 
127 


where C2(G) and T(R) are Casimir operators for the adjoint representation G of the 
color group and for the representation R of the fermions respectively. In the case of 
SU(\) color With f flavours, they are given by 


; (14) 


C2(G) = Cabe Cabe = N , (15) 


N? — 1 abc 


TCR) Sap = Tr(¢4t”) =$ fap - (16) 
For the validity of eq. (13) it is necessary that Q3, in eq. (10) be chosen large 
enough to make a(0)/r sufficiently small for a lowest formula to apply. Finally 
ANS are a set of constants independent of œ which we shall study in the following. 
It is immediately seen that, for each n, eq. (12) is the solution of the differential 
equation 
dM S(t) _ = 
dt 
with assigned initial value MNS = MNS(0) at t= 0. In turn, for any n (sufficiently 


large) the whole set of eqs. (17) is equivalent to the following master equation for 
the densities: 


ae t) == 2 jg So, 1) PÈ), (18) 


DA ANS MNSGA), (17) 


provides that 
1 
fa z”! pz) = ANS | (19) 
0 


This is easily been by multiplying both sides of eq. (18) by x”~!, integrating from 
zero to one and interchanging the order of integration on the right-hand side. Eq. 
(18) is our starting point since it allows a simple physical interpretation of the func- 
tion P(z) which we shall use in the following for its direct computation (and of 
other similar generating functions). 
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It is convenient to rewrite eq. (18) in the form 


1 1 
qNS(x, t) + dqNS (x, t= fay a &(zy — x) qS, Dic —1)+ z5 P(z) ar . 
(20) 


The meaning of this equation is clear. Given a quark with momentum y there is a 
chance that is radiates a gluon, thus reducing its energy from y to x. The quantity 


Paa + dPaq = (Z — 1) += PG) dt (21) 


is the probability density of finding, inside a quark, another quark with fraction z 
of the parent momentum. The change with t of this probability produces the varia- 
tion of the quark distribution function. Thus A(z) a/27 is the variation per unit ¢ 
at order & of the probability density of finding inside a quark another quark with 
fraction z of the parent momentum. 

We now drop the restriction to one flavour and to non-singlet densities. In parton 
language a singlet density is in general a combination of the sum of all quark and 
antiquark densities and of the gluon density inside the proton. We therefore intro- 
duce G(x, t) as the density of gluons (summed over colors) inside the proton in the 
Po frame. We can now directly write down the integro-differential equations that 
describe the Q? dependence in the general case. They are 


a = jel 2 Jy, Pata! (=| + GQ, t) Pio(*)| f (22) 
sa, 2- x0 f a] x q'(y, t) Pegi (=) + Gy, t) Pec (*)| : (23) 


Here the indices i and j run over quarks and antiquarks of all flavours. The number 
of quarks as seen by the current changes by two mechanisms: a quark originally at 
higher energy may loose momentum by radiating a gluon, or, a gluon inside the pro- 
ton may produce a quark-antiquark pair. Similarly the number of gluons changes 
because a quark may radiate a gluon or because a gluon may split into a quark-anti- 
quark pair or into two gluons. This last possibility is typical of non-Abelian gauge 
theories where a three gluon vertex exists to order g (while the four gluon vertex 
is of order g°). 

Some properties of the functions P(z) appearing in eqs. (22), (23) are immediately 
derived from the fact that color and flavour commute. First Pyig/ is diagonal in 
quark indices because a gluon is emitted without flavour exchange, 


Pgigi = 5 i Paa ; (24) 


Moreover, when we neglect all masses, the probability of emitting a gluon is the 
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same for all flavours, 
Paa! = Pog (independent of i) . (25) 


Finally a gluon creates a massless quark-antiquark pair with equal probability for 
all flavours. Thus, 


Paic = Pac (independent of i) . (26) 
Therefore we can rewrite eqs. (22), (23) in the simpler form 
dq‘(x, t N} i 
aie ) = i a O, £) Py E) GU, 1) Pac(*)| (27) 
dG(x, t) a(t) pdy a x x 
T = an f: y [£ q'(y, t) Peal] + GO, t) Pool) (28) 


By summing eq. (27) over i = 1, ..., 2f, we obtain 


dE? q(x, t) O fD [2 d0, 1 Pale )+ G0 t) 2f Pao ~)]. (29) 


dt 20 be i=1 
It is precisely the matrix 
1 Pa) SPa) | |ANS 4R) 
dz z” ~? = : (30) 
0 Pea) Pac) AŞI ARS 


which gives the logarithmic exponents for each n as given in the literature [4,5]. For 
each value n this matrix must be diagonalized in order to find the eigenvectors and 
the eigenvalues of the Q* evolution equations. 

By subtracting the derivatives for a quark and an antiquark (or for two quarks) 
we immediately recover the equation for the non-singlet case, eq. (18), 


$ ies D — ae, D) 


= a fe BOF Paa(= ) | (31) 


We thus reproduce the result that Pg, is the same function for singlet quarks and 
for all types of non-singlet quarks. 

The non-diagonal functions (a/27)Pgq(z) and (a/27)P,g(z) can be directly inter- 
preted as probability densities (a/27)Pgq(z) is the probability density per unit ¢ at 
order & of finding a gluon inside a quark (or an antiquark) with fraction z of the 
longitudinal momentum of the parent quark. (a/27)Pc(z) is the probability density 
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per unit ¢ at order a of finding a quark (or an antiquark) inside a gluon with fraction 
z of the longitudinal momentum of the parent gluon. For the diagonal functions 
(a/2m)Pgq(z) and (a/2m)Pgg(z) the actual probabilities also involve delta function 
singularities at z = 1. The probability densities are given in these cases by eq. (21) 
and by 


Poa +d Pog = 6(z — 1) + 5 Poc(z) dt , (32) 


Therefore it is only at z <1 that (a/21)Pgq(z) and (a/27)Poc(z) are probability 
densities (for example they must be positive defirtite). Consider for example 
(a/27)P,q(z). Since the total number of quarks minus antiquarks is conserved, the 
probability of finding a quark in a quark, integrated over all values of z, must add up 
to one. It follows that the integrated correction of order œ must be zero 


1 
faz Paa) =0, (33) 
0 


which is a condition on ANS, This is the well known result that charges (i.e. current 
algebra sum rules) are protected against Q? corrections [4,5]. We see that the pres- 
ence of ô function singularities destroys the positive definiteness of integrals of Paq 
and Pgg. In particular we shall see that the values of ANS are all negative for n > 1. 

Momentum conservation in the vertices imposes further constraints on the P 
functions. At z < 1 we have 


Pag (2) = Pgg(l =z); 
Pac (z) = Pag — 2) (z<1), (34) 


Pgc(2) = Poc! — 2). 


The above equations arise because when the quark radiates, it splits into a quark 
with fraction z of its momentum plus a gluon with fraction (1 — z) etc. However at 
z = 1 the 6 function singularities upset these relations. But it must remain true that 


1 
faz 2[Paq(2) + Poa) =0, 
7 (35) 
faz zif Pyc@) + Poc) = 0, 
0 


which guarantee that the total momentum of the proton (i.e. of all partons) is un- 
changed, 


ae ae 
a f dx x[ 2 qi(x, t) + G(x, »|=0 (36) 
0 Iz 
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as can be seen from eqs. (28) and (29). In the following we shall use eqs. (33) and 
(35) to fix the 6 function singularities at z = 1 of Pyg and Peg. 

We can make the symmetry relations eqs. (34) hold for all values of z if we 
treat the points z = 0 and z = 1 symmetrically i.e. if we artificially add (for the sake 
of this argument only) a 5 singularity at z = 0 when the symmetry requires it. Note 
that moments with n > 1 are unaffected by a 6(z) term. The interest of this obser- 
vation is that we have, for example, 


1 1 
fa zh". (2) = ANS = fa z"—! Pgg(1 —z) 
0 0 


1 n—i 
r _ n-i £ er ae ec 
feo z)"—! Poq(z) 2 ¢ 1) ("" ) 483, . (37) 


Through eq. (33) the previous sum rule would imply for n = 1 that Aga = 0. This 
equality is not true but is the only consequence of having artificially added a 5(z) 
term to Pgq. We can however use it to eliminate asa from eq. (37) and write a sum 
rule for the logarithmic exponents with n > 1 which reads 


n—l 
n-1 
ANS = 2 (-1)* aSa. (38) 
k=1 k 


This sum rule must be (and is) satisfied by the calculated logarithmic exponents 
and offers a general test of the above physical interpretation. The origin of eq. (38) 
in the light cone language is not easy to disentangle. Another similar sum rule is ob- 
tained for the combination 


By = 2f ATS + AGS (39) 


by using the symmetry relation eqs. (34) and the fact that B3 = 0 as stated in eqs. 
(35). This sum rule reads 


n-1 
n—i 
B, = 2 (-1)* ( ) Bk+1 > 
k=2 k 
or equivalently, 


n- 
-u ni 
(+01, = D ei (™ ) a. (40) 

j=3 j-1 

This relation implies that the B,,’s for even n are determined by the B,,’s for odd 
n. 

We conclude this section by observing that master equations of the form of eqs. 

(27), (28) are valid in any renormalizable theory at lowest order in perturbation 
theory. Therefore they are reliable provided the coupling constant is sufficiently 
small. For example they are valid in quantum electrodynamics of electrons and 
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photons. In this case Pag > Pyy ~ 5(z — 1) because of the absence of a 3y coup- 
ling. The replacement of a by a(t) in QCD guarantees that for sufficiently large t 
the master equations will eventually become reliable. These master equations can 
provide a useful insight in connection with problems that cannot be approached by 
light-cone expansion. 


4. Calculation of logarithmic exponents. Spin averaged case 


In this section we show that the functions A(z) introduced in the previous section 
can be directly computed from the simple knowledge of the basic vertices of QCD. 
The method used is an extension of the von Weizsacker-Williams result in quantum 
electrodynamics [8]. In that case the equivalent number of photons inside an elec- 
tron with fraction z of the electron momentum is evaluated to order œ and contains 
a factor of In E/m., which plays the same role as ¢ = In Q7/Q2 in our case. 

We first compute the functions P(z) at z < 1 while we shali deal, at the end, with 
ô function singularities at z = 1. We want to evaluate the probability of finding a 
particle B inside a particle A with fraction z of the longitudinal momentum of A in 
the p, frame to lowest order in a: 


dP pa(z) dz = = Pya(2) dz de. (41) 


Let C be the third particle in the bare vertex where A and B appear. We can 
identify the above probability by comparing the cross sections for the two processes 
in fig. 1, where D is a given particle and f any final state. We define the general S- 
matrix element as 


Sy — öy = 2mi8 (Bj ~ E) My TT QED} , (42) 


where the index k runs over all external particles. Although by no means necessary, 
we find it particularly useful to phrase our calculation in terms of the “old” pertur- 
bation theory which is best suited for a discussion of leading terms in the p,, frame. 
The contribution to Mj; in eq. (42) of a given intermediate state B to the process in 


Fig. 1. (a) Contribution of the B intermediate state to the process A+ D >C + f. (b) The pro- 
cess B+ Df. 
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fig. 1a can be written as 


Vaspic VB+D>f 
® TOE Ep t Ec — Ea)’ 


where V,; is the invariant matrix element of the interactions (with the factors (2E,) 
removed). Similarly in the case of fig. 1b we have 


Ma+p-c+t = (43) 


-1/2 


Mpip—t=& VB+D>f . (44) 


The contribution of a given intermediate state B to do, and doy, is therefore given 
by 


2 
g? WVa+specl?|Vp+p—el 


doa = 8B Ep QE) (Ep + Ec - Ea)? 
254 ake dpe 
X (2m) ô (ka tkp- ke ~ ke) OR Qn QED’ (45) 
do, = BE we IVa+p—+el7(20)*54(kptkp— koll aati. (46) 
The two processes are related by 
dog = dP ga (z) dz doy . (47) 
From eq. (47), by the comparison of eqs. (45) and (46) we obtain 
2 2 3 
tO anh =F OE Mba be Ex) COED e 
where all masses being neglected: 
ka =(P; P, 0), 
2 
kg = (2p + 2 SP’ ;2P, pı). 
: pi 
o (0 -2P+ Ps -2P -p1). (49) 
We thus have 
2 2. pP? 
(2EB) (Ep + Ec — Ea) ere (50) 
deke _ dzdp? 61) 


mP QE) 16r?(1 — z) ` 
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The result is 


a 21 —z) > Va-+pecl? 
2 


dinp’, (52) 
2r 2 spins Pi Pi 


dP ga(Z) = 
where a sum over the spins of B and C and an average over the spin of A is indicated 
(if the case). 

In the case of interst |V]? vanishes linearly in p? so that the ratio |V|? /p? in eq. 
(52) is finite at p? = 0. Moreover for a virtual mass —Q? for particle D, the integral 
in p? has an upper limit of order Q?, so that, at the leading logarithmic approxima- 
tion, d In p? can be directly interpreted as dt. We thus finally obtain, by comparing 
eq. (52) with eq. (41), 


a 2 
Poa@)=$ 201-2) D Warpecl gxi), (53) 
spins Pi 


an expression which only depends on the vertex ABC. In particular when the spin 
sum is symmetric we clearly have 


Poa(z) = Pga (1 — 2) (z<1) (54) 


in agreement with eqs. (34). 
Let us now specialize to the quark gluon vertex in fig. 2 in order to evaluate 
Pgq(z). In this case, 


D asoa = $Ca(R) Toray) Hete, (55) 


spins ° 
where the factor of 5 derives from the average over the initial quark spin, and C2(R) 
defined by 
N?_] 
2N 
arises from the sum and average over the final and initial states in color space. Care 


must be taken so that only physical transverse gluon states are included in the sum, 
and we therefore write 


C,(R) = iD “p = (56) 


o kiki 
Drete >i TBB  (ij=1,2,3). (57) 
pol k5 


One obtains 


Fig. 2. The quark gluon vertex which determines Pga and Pgq: The form of the vertex is 
iggcy'ttq ABI, with Trt? = 359°. 
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Fig. 3. The annihilation vertex of a gluon into a quark-antiquark pair which fixes Pgg.- 


2p? 1+(1-z)? 


Va+c+ql? = (R). 58 
2| C = 2(R) (58) 
We can thus state the result from eqs. (53), (58), 
1+(1 —z)* 
Peale) = CR) 7 E— (59) 


which holds at all z, since we are dealing with a non-diagonal density. From the last 
equation, by using the symmetry relations eqs. (34) we also obtain 


2 
PaO ON <1) (60) 


We now calculate Pyg from the vertex in fig. 3. Since Pqg is proportional to the 
probability density of finding inside a gluon (averaged over colors) a quark (or an 
antiquark) of given flavour and of any color, in this case the sum and average in 
color space simply bring in a factor of A. We thus have 


D Weoqeal? =} Hc ruta) } Dette, (61) 
spins - pol 
which, recalling eq. (57), gives 
l-z 
D osar? =a no ) 


spins l-z 


From eq. (53) we finally obtain 
Pac(2) = 32? +0- 2)°). (62) 
The symmetry under the change of z into (1 — z) is expected because Pyg(z) 


= Pag(1 — z) = Pag( — 2), by eqs. (26) and (34). 
We now turn to the three-gluon vertex in fig. 4 which determines Pgg(z). The 


A by 


ap 


cÀ 
c 


Fig. 4. The three-gluon vertex relevant to determine Pgg. It is equal to 


—igcane(Sypka + kB)a— Surtke +ka)v t+ grvkkc — kp)ul - 
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amplitude is given by 


Vo+c+G = —Cabe{—[ka + kc) e$ 1 (€4 €6°) 


nE [kc - —kpeA] (ec ep?) + [Kat kg)ec'] (GA en yh (63) 
It is only a matter of algebra to derive the result * (recall eq. (15)) 


spins 


eee Éz 
L lVasexel? =4C,(G) a A “tig |. (64) 


Through eq. (53) this equation leads to 


Z 
1 = 


Note the obvious symmetry for z into (1 — z) in agreement with eq. (34). 

We now complete the determination of Pg,(z) and Pgg(z) by fixing their be- 
haviour at z = 1. Note that all moments of these functions as given by eqs. (60) and 
(65) would be divergent at z = 1. We therefore start by regularizing the factor 
(1 — z)~! by reinterpreting it as a distribution (1 — z);! defined in the following 
way: 


Pools) =20(6)| 2 + ~ + (1 a| (2<1). (65) 


rafo _ A IO- E, d 
laa eae ea (66) 


with f(z) being any test function which is sufficiently regular at the end points. In 
particular we have 


i 1 
d — z 3 6 
Jen 0 (67) 


We then add, to Pgg(z) and Pgg(z), a 6(z — 1) function with the coefficient 
determined by the constraints in eqs. (33) and (35). We thus find 


2 
Peal) = CR) | A= +2o@—n), (68) 
Poc (z) =2C2(G) i A: + a +2(1—z)+ (4 = Tra 6(z — D| i 
(69) 


It is now straightforward to evaluate the moments of the P(z) functions, which 
give the set of constants A, according to eq. (30). We first evaluate the moments of 


(1—2z);'; 


* We shall study this amplitude in all details when dealing with the polarized case in sect. 5. 
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1 gl 1 zy _ 5N (-1)/ n-i 
Nea (ear Seas l-z -2 ; ( ’) 


j=l 1 i 


n=l 


OS EA (70) 


j= 


_ 
~, 


where the last equality can be easily proved by induction. We finally obtain 


1 1 a 

dg z"~! = ANS =C (R) |-4+——. -22, >], 71 
fee As A ear coer A 
- ki _— 4G 2tntn? 
- dz z"~! Pg (z) = AS4 = CAR) GF f (72) 

| 2+ntn? 

n-1 = G- Sa ante eee 

ue Pac) = RAs ZR) ein #2)’ Ca 
1 
fu z"-1 Pool) = AGS 
0 

i I 2 2 1 2 TR 

SCG) |- N n(n —1) y (nt1)n+2) “j=2j 3 ao) oe: 


This set of logarithmic exponents, taking eq. (12) into account, is seen to coin- 
cide with the results of refs. [4,5]. 


5. Calculation of logarithmic exponents. Spin-dependent case 


In this section we shall consider the Q? dependence of quark densities with 
given helicity which are relevant for scaling breaking effects in deep inelastic scatter- 
ing on polarized targets. 

The experience accumulated over the previous sections allows us to directly write 
down the general master equations which are 


£ gi, t)= © fe [éo DPaae(=] +4 0, DPa,a (>) 


sear orale) ooon], 
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£ gl (xt) = 2 Fao. t) Pg_ a-(5)* 40,1) Paal) 
+60, o,_c,(=) PEO Dro (5) i 
yL 


F a(t) 1 at x 2 
SE i 2 i x 
= Gr) =o f |: ES ab È gib; d Po,a(>) 


+ GQ, t) Po,¢,(=)+6- Q, ) Pec 5)) 


soan JOD dodro af) oora) 
+640, 0) Po_c,(*)+6-0. Poo. ~|I. (75) 


The previous set of equations can immediately be simplified by observing that 
parity conservation in QCD implies the relations: 


Pa B,C) =Pa_s,(2) (76) 


for any A and B. These relations in turn imply, (when used in eqs. (75) that the 
sums q} + qi =q' and G, + G_ =G and the differences 


Ad =a -q', (17) 
AG=Ĝ; -G ; (18) 


evolve separately. For the sums we recover eqs. (27) and (28) with the obvious 
identification that 


Par = Pa B, + Pa_s, - (79) 
It is convenient to define 
APAs =Pa,p, — Pa_s, - (80) 


We then derive, from ja (75), the master equations for the differences 


Š Aai(x, D= a Z fagio, #) APao(=)+ AGO, D APyo(=)] 


$ AG(x, t) = a0 i ap Aq'(y, t) APoq E) + AG(y, t) arec(=)| - (81) 
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In the simplest case of non-singlet quark densities, eqs. (81) reduce to 


o 


5 Aas. D fF aSo, 2) APaf >). (82) 


For example this allow us to conclude from eq. (8) that the non-singlet part of 
the structure function G, is an eigenvector of the Q? dependence. Actually we can 
immediately obtain a much stronger result in that, when masses are neglected, the 
vector quark-gluon coupling in fig. 2 is helicity conserving. Therefore we can direct- 
ly conclude that 


Py_a,(2) =0, 


Paya, ©) = Peal) = Pa = CAR) 1 
We thus obtain that the moments of qNS(x, t) and of AqN5(x, t) (or equivalently 
of the non-singlet parts of F, and G,) evolve in Q? with the same logarithmic expo- 
nents. 
The same quark-gluon vertex of fig. 2 also determines Pg, g,(z) and Pg_q,(z). 
Given the gluon momentum in the form (see eqs. (49)) 


pa 38e- 1]; (83) 


zP 


for the corresponding polarization four-vector for positive and negative helicity we 
can approximately set 


Px + ip ; 
e+ = (0; — 2 viivi), 85) 


where only the leading term for each component was kept, which is sufficient for 
our purposes. For the invariant vertex Vy,+G,+¢ (see eq. (55) for comparison; a 
sum over the final quark spin is understood) one has . 


p? 
= (r + oP ; ZP, Px, Ps) (84) 


l+y 
Vo,>Gaql? = C2(R) Tr(Hequat 5 s) ej"e; 


ae 


ES Lie- 2z +z?) + z(2 — z)] . 


By using eq. (53) we thus obtain 


Pory D =R, (86) 


2 
Po_ai(2)= CR) =, (87) 
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and consequently, 


_ _ pje 
APgaq(z) = C2(R) ae i (88) 


while the sum of the two reproduces eq. (59) for Pgg(z) in agreement with eq. (79). 
We now evaluate P,, g,(z) and Pg_g_,(z) from the annihilation vertex in fig. 3. 
A sum over the final antiquark spin being understood, the invariant vertex gives 


Voy-asal” = 5 Te(torston ! =1s) A 
pi 2 2 2 2 
"aa Erao -0-9)1. (89) 


(Note that €„ now refers to the initial gluon so that it is given by eq. (85) with 
pı = 0). From eqs. (53) and (89) we obtain 


Pua) = 427, (90) 
Pac) =3(1 -2 , (91) 
APD = 3 [27 -(—z)7], (92) 


while the sum reproduces Pag in eq. (62). Last but not least, we must now evaluate 
Poca, and Pg_g, from the three gluon vertex in fig. 4. The invariant amplitude 
was written down in eq. (63). When the impulses of the gluons are specified as in 
eq. (49) the polarization vectors are 


ELA = J1(0; 0, 1,i) ? 


Py, ti 

en = VA(0;— x= Py, xi), (93) 
ti 

Esc -io Ee, 1, +i) 2 (93) 


It is more expedient to calculate the four helicity amplitudes first and to take the 
square modulus and sum over the polarizations of the gluon C at the end. A simple 
calculation leads to the following results for Vg,+G4G4' 


, 1 1 
Vi++ = V2 Capc(Px — ipy) (; + h} 


1 
V}-+ = -v2 CabePx + ipy) (;- 1) > 
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1 
Vit = —V2 Cane(Px + ipy) (=; = 1) > 
V,-_=0. (94) 
By using eq. (53), a simple algebra leads to 
1 
Po.64@) = C,(G) (1 + z4) (- + ==). (z<1), (95) 
a-z)? 
Po aO CG) E. 06) 


The behaviour at z = 1 of the diagonal density Pg, g,(z) is immediately obtained 
from that of Pgg(z) in eq. (69) by recalling eq. (79) (which is verified by eqs. (95) 
and (96) at z < 1). We thus obtain 


l jes TN Joe - Don 


PaO CHG) [+29(2 +7 


Z)4 z 


aroa Olalia C= -iR e-o] 
(98) 


Collecting our results eqs. (83), (88), (92) and (98) we are now ready to compute 
the matrix of the logarithmic exponents 


1 APyq(z)  2fAPyc(z)| | ANS  47(R)A2° 
f dzz"~! y =| ; (99) 
ò APcq(Z) APool) 4Sa ASS 
which are 


JNS — 4NS 
An =A, ? 


-1 
ne’ 


47(R) AGS = 27(R) 


n-i 
R a 2D 4] (100) 


ll 
ASS = €,(G E ~ 
CXS) 3GG@ n nè “jaf 


Taking eq. (12) into account, this set of values is in agreement with the result of 
refs. [6,7], obtained by the operator formalism. 
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6. Concluding remarks 


We briefly comment on the reasons for the success of this simple method of 
calculating the logarithmic exponents and the limitations of the present approach. 
The essential point is the following. The exact renormalization group equations are 
given by 


a a x 
[2 -1o 2 -m| ma )=0. (101) 


We refer to the simplest case of a multiplicatively renormalizable operator, whose 
coefficient function is proportional to the moment M,, of some structure function. 
At lowest order in a one has 


Q 
n(@) ~ = An - (102) 
n 


Eqs. (12) and (13) show the solution of eq. (101) to lowest order in œ. From eq. 
(17) we see that the solution of eq. (101) also satisfies the simpler equation to low- 
est order in a, 

d _| daft) os 
È = ralo()| Mn -|$ r. O 4,] Mn =0. (103) 

Thus, to lowest order the only effect of the 8 function, which embodies all the 
intricacies of the coupling constant renormalization, is to teplace a by the running 
coupling constant in the expression of y,,. Therefore the lowest order master equa- 
tions through this simple replacement become the solutions of this more sophisti- 
cated problem. In higher orders the effects of (a) cannot be taken into account 
in this trivial way and the evaluation of the functions A(z) would be more com- 
plicated. 

Another feature which makes the calculations particularly simple is the direct 
and transparent connection of the relevant bilocal operators in terms of parton 
densities of given helicity in the p. frame. It would be interesting to investigate 
whether a similar approach can be applied to the computation of the logarithmic 
exponents of more complicated operators, as those, for example, which are relevant 
to the structure function G2. 
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